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DIFFUSION OF HEAT FROM A LINE SOURCE IN ISOTROPIC TURBULENCE ^ 


By Mahinder S. Uberoi and Stanley Corrsin 


SUMMARY 

x\n experimental and analytical study has been made of some 
features of the turbulent heat diffusion behind a line heated wire 
stretched perpendicular to a flowing isotropic turbulence. The 
mean temperature distributions have been measured with sys- 
tematic variations in wind speed, size of turbulence-producing 
grid, and downstream location of heat source. The nature of 
the temperature fluctuation field has been studied. 

A comparison of Lagrangian and Eulerian analyses for dif- 
fusion in a nondecaying turbulence yields an expression for 
turbulent-heat-transfer coeflicient in terms of turbulence velocity 
and a Lagrangian ^ ^ scale. 

The ratio of Eulerian to Lagrangian microscale has been de- 
termined theoretically by generalization of a result of Heisenberg 
and, with arbitrary constants taken from independent sources, 
shows rough agreement with experimental results. 

A convenient for?n has been deduced for the criterion of inter- 
changeability of instantaneous space and time derivatives in a 
flowing turbulence. 

INTRODUCTION 

One of the most striking aspects of turbulent motion in 
fluids is its dispersive property. This ^ ^convective diffusion/^ 
illustrated by the general statistical tendency of (noncon- 
tiguous) fluid elements to get farther apart with increasing 
time, was probably first observed long before the era of 
analytical fluid mechanics. An analytical start on this 
problem was not made, however, until the now-classic work 
by Taylor in 1921 on difl‘usion by continuous movements 
(reference 1). Not only did this paper lay a groundwork for 
the study of turbulent diffusion but it also represented a 
forward step in the ideas essential to development of a gen- 
eral statistical tlieory of turbulence, a field which had scarcely 
progressed since Reynolds’ original formulation of the equa- 
tions of motion for a flow in which mean and fluctuating 
parts could be distinguished. 

The diffusive action of a turbulent flow may manifest itself 
in various ways, depending upon the initial and/or boundary 
conditions and upon the interests of the observer. The fol- 
lowing possible measures of the diffusive powers are neither 
exhaustive nor mutually independent: 

(1) The average rate of dispersion of particles from a fixed 

source 

(2) The average rate of increase of spacing between differ- 

ent particles 


^ Supersedes NACA TN 2710, “Difiusion of Heat From a Line Source in Isotropic Turbii 
lence” by Mahinder S. Uberoi and Stanley Corrsin, 1952. 


(3) The average rate of transport of particle concentration 

under a given mean concentration gradient 

(4) The average rate of increase of the length of a fluid line 

(5) The average I’ate of increase of the area of a fluid sur- 

face 

The word '^particles” means simply indelibly tagged fluid 
elements, much smaller than the smallest length associated 
with the turbulence. 

The present report is concerned primarily with measure 
(1). The measurements have all been made in the thermal 
wake of a long thin heated wire mounted perpendicular to 
an isotropic turbulent air flow and producing no turbulent 
wake. Here the tagging is tliermal, and the degree of in- 
delibility (negligibility of molecular diffusion) is one of the 
matters to be investigated. 

The difl’iisive property (for scalars) of a turbulent flow is 
apparently a secondary characteristic at least in the sense 
that it need not explicitly enter the dynamical problem. The 
diffusion may be regarded as a kinematic phenomenon, to be 
deduced from the dynamical solution to the problem if and 
when the latter is obtained. Thus the objective of research 
on turbulent diffusion may be to seek a connection between 
the diffusive and the dynamical statistical variables, even 
before the complete dynamical theory is available. 

Measure (3) is usually termed the ‘turbulent transport” 
or ^ transfer” problem. Although of extensive practical im- 
portance, it has not yet been subjected to genuine theoretical 
study. 

Most of thesemiempirical “theories” of turbulent transport, 
for both scalar and vector properties, employ an Eulerian 
formulation of the basic equations, and up to now they have 
been unable to relate the turbulent transport correlation to 
other statistical functions describing the flow. Taylor 
(reference 1) showed that in the simple case of a homogeneous 
field of isotropic turbulence, and even in a decaying isotropic 
turbulence (reference 2), a Lagrangian formulation of the 
transport (i. e., diffusion; the terms will be used interchange- 
ably) problem leads to some important results.^ 

Up to the present time little theoretical or experimental 
work has been done to find relations, if any, between the 
Lagrangian statistical measures of a turbulent field with its 
Eulerian statistical measures. Since turbulence dynamics 
seems best handled in the latter terms and turbulent diffusion 
in the former, it is evident that such a connection is impor- 

2 In his tensorial generalization of Taylor’s work on case (1), Batchelor (reference 3) has 
chosen to call this an “Eulerian” analysis, describing only case (2) as “Lagrangian.” In 
keeping with previously accepted nomenclature, both cases are Lagrangian, (1) involving a 
single particle and (2) involving a pair. In fact, case (2) might be termed a mixed (Eulerian 
and Lagrangian) problem. 
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tant. Hence, one of the purposes of the present experiments 
has been to compare the magnitudes of some of these quanti- 
ties under variations in the turbulent field. For example, 
the postulates of Taylor and Heisenberg on a relation between 
Lagrangian and Eulerian microscales can be examined and, 
in corrected form, compared with experiment. 

The turbulent dift'usion from a fixed line source can be set 
up analytically as an ordinary (Eulerian) “heat-transfer” 
problem, permitting a start to be made in relating measures 
(1) and (3) of the diffusive power of a turbulent flow, under 
certain simplifying assumptions. 

Measures (4) and (especially) (5) may well be classed as 
characteristic of the “turbulent-mixing” problem rather than 
diffusion in the common connotation. 

Experimental work on diffusion from a fixed local source 
in a turbulent flow has been meager. In isotropic turbulence, 
there have been the measurements of Schubauer (reference 
4), Simmons (reported by Taylor in reference 2), Dupuis 
(reported by Kampe de Feriet in reference 5), Frenkiel (refer- 
ence 6), and Collis (reference 7). Of these, only the data of 
Simmons and Collis are extensive enough to permit confident 
computation of the Lagrangian correlation function. In 
turbulent shear flow, Skramstad and Schubauer (reference 8), 
Dryden (reference 9), and the present authors (reference 10) 
have measured distributions close to a source; Kalinske and 
Pien (reference 11) and Van Driest (reference 12) have made 
measurements somewhat farther downstream. 

None of these studies was repeated with a systematic 
variation of the properties of the turbulence. In spite of 
the poor precision inherent in this type of measurement, it 
was hoped that such an approach would at least show up 
some general trends in the relations between Eulerian and 
Lagrangian variables. 

This investigation has been conducted at the Aeronautics 
Department of the Johns Hopkins University under the 
sponsorship and with the financial assistance of the National 
Advisory Committee for Aeronautics. The authors would 
like to acknowledge the assistance of Messrs. Alan Kistler, 
George Stierhoff, and Allen Gates and Miss Patricia O’Brien, 
as well as the helpful criticism of Dr. Francis II. Clauser and 
Dr. C. C. Lin. 

SYMBOLS 


c root-mean-square molecular velocity 

Cp specific heat per unit volume at constant pres- 

sure 

f{r) Eulerian velocity correlation coefficient (nota- 

tion of Von Karman and Howarth) 


H = 2pc 



//*- 


2pc 


j 

k 

k 7 ’ 

L 


II 

n 00 

Jo 

width of rectangular heat pulse 
thermal conductivity 
turbulent-heat -transfer coeflScient 


Eulerian scale 




Lagrangian scale (^Ll ■j; Rv{r)dT^ 
Lagrangian scale for nondecaying and decay- 
ing turbulence -j; Hr,{ri)dT^ 

mixing length 
grid mesh size 

= 7jRr,['n)dr] 

P{d^) probability density of fluctuating temperature 


Ll 


I 

M 


V 

Rl 


static pressure 

turbulence Revnolds number based on Euler- 


ian scale 




P = P (f I "I v{t)v{t r) 


Rv{t) 

Rih) 

p 


Rx 


Lagrangian correlation coefficient for nonde- 
caying turbulence 

Lagrangian correlation as a function of rj for 
nondecaving and decaving turbulence 


turbulence Revnolds number based on Euler- 


ian microscale 




r scalar distance between two points 

--(U/[(iy]” 

s average on-center spacing of pulses 




U 

u 


u' =\U^ 


V 


time 

mean velocity in ^-direction 
instantaneous velocity fluctuation in a:-direc- 
tion 

mean velocity in ^-direction 
instantaneous velocity fluctuation in y- 
direetion 


w 

X 

X 

To 


instantaneous velocity fluctuation in z- 
direction 

distance traveled in j:!-direction by a fluid 
particle 

distance downstream from grid 
location of heating wire 


\x = x — x 

Y 

U' 

Ym 

V 


distance traveled in y-direction by a fluid 
particle 

root -mean-square displacement of a fluid 
particle in y-direction 

root-mean-square displacement of a molecule 
distance in direction of measured diffusion 
distance in direction of heating wire 
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n(^oiO= f 

O instantaneous temperature (measured above 

ambient room temperature) 


0' = -yld'^ 

K 

A 

X 


X. 

X. 


V 


mean temperature 

maximum mean temperature at a cross sec- 
tion, a function of Aj:- 
instantaneous temperature difference 


temperature difference of rectangular heat 
pulse 

dimensionless empirical constant 
mean free path of a molecule 
Eiderian microscale of turl)ulence 

Laerangian microscale of turbulence for non- 

/ / 2 

decaying turbulence = "(6) ) ) 


J^agrangian microscale of turbulence for 
nondecaving and decaying turbidence 

kinematic viscosity 


p 

T 

(-) 

Subscripts: 

max 

min 


density 

time difference 

mean value or ensemble average 

maximum 

minimum 


EQUIPMENT AND PROCEDURE 

AERODYNAMIC EQUIPMENT 

The wind tunnel (fig. 1) is an open-return NFL type tunnel 
with a 2- by 2-foot working section and a^ree-stream turbu- 
lence level of v'IU=0M percent and u'IU=^M percent at 
a mean velocity of 26 feet per second. The turbulence- 
producing grids were as follows: 


Designation 

Type 

Mesh size 
(in.) 

Rod 

diam- 

eter 

(in.) 

Solidity 

1-in. grid 

^-in. grid — 
H-in. grid-... 

Biplane, wood.. 
Woven, metal . 
Woven, metal. . 

1.00 

«. 50 by 0. 53 
.25 

0. 25 
. 120 
.060 

0. 437 
.41 
.42 


« The 0.50-ln. mesh was set in direction of measured diffusion. 



5-hp blower 


3 *-] 


-7'- 


Grid 

Square to round transition sections ' 



Over-oil length « 42^ ft 
Figure 1. — Sketch of open-return wind tunnel. 


Cheesecloth screen 

All other screens - #24 wire mesh 


They were mounted in turn at the upstream end of the work- 
ing section. 

The heat source was an 0.008-inch-diameter platinum wire 
stretched vertically across the tunnel at various distances 
from the grid. It was heated by direct current to tempera- 
tures between 500° and 700° C, with the latter figure only 
at the highest operating velocit}" of 38.0 feet per second. 
The wire Reynolds numbers at this condition and at the two 
other velocities were as follows: 


Mean 

speed 

(fps) 

Wire 

tempera- 

ture 

(“C) 

Reynolds number 
based on— 

Air tem- 
perature 

Wire 

tem- 

perature 

8.5 

500 

35 

6.4 

25.6 

500 

105 

19 

38.0 

700 

157 

19 


A prelimiiiaiy investigation was made without grids to 
insure that these operating conditions did not generate a 
vortex street downstream of the heated wire. 

With the grids in place, the mean momentum wake became 
practically undetectable with total-head tube and manometer 
at distances greater than 1 or 2 inches downstream of the 
heating wire. 

MEASURING EQUIPMENT AND PROCEDURES 

The mean-temperature distributions were measured with a 
Chromel-Alumel thermocouple and a Leeds and Northrup 
type K-2 potentiometer. The cold junction was kept out- 
side of the wind tunnel. 

The shadowgraph technique was used to photograph the 
laminar thermal wake close to the source with no grid in the 
wind tunnel. This information was applied to the problem 
of “correcting’’ the thermal wake in turbulent flow for the 
eflects of molecular diffusion and of finite source size. Figure 
2 (a) is a shadowgraph of the wire wake with no grid in the 
tunnel; figure 2 (b) is a typical time exposure with grid- 
produced turbulence. A resistance-thermometer traverse of 
the laminar wake in a flow of very small turbulence showed 
that the temperature profile had already become very nearly 
Gaussian at a distance of 1 inch (125 wire diameters) down- 
stream. The white lines on the sides of the dark wake 
shadow in figure 2 (a) correspond to the minimums in the 
second derivative of the density profile. For small tempera- 
ture differences these coincide with the maximums in the 
second derivative of the temperature profile. Although the 
temperature differences are not small in the immediate 
vicinity of the wire, this condition is reasonably well satisfied 
at relatively small values of Ax as evidenced by the parabolic 
spread of this laminar wake. 

The standard deviation of the wake could then be com- 
puted from the spacing of these two bright lines, with the 
assumption of a Gaussian distribution. Since the closest 
points of traverse in the turbulent cases were K inch (63 wire 
diameters) from the heat source, this was probably a reason- 
able assumption. 

As pointed out by Taylor (reference 2), the molecular and 
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(a) No grid in tunnel. 

(b) Grid-produced turbulence. 

Figure 2. — Shadowgraph time exposure of wire wake. 

the turbulent diffusive phenomena are statistically inde- 
pendent, so that the squares of the standard deviations due 
to these two effects are additive. Hence the wake spread 
due to turbulence alone, from a true line source, was obtained 
by subtracting the square of the standard deviation of the 
laminar wake (computed from the shadowgraph) from the 
square of the standard deviation of the total wake (from 
thermocouple traverses with grid turbulence present) at all 
stations. This difference was the square of the standard 
deviation of the desired phenomena. All wake-spread data 
presented in the next section have been corrected in this 
fashion. 

Parenthetically, it should be remarked that the laminar 


wake in figure 2 (a) spreads parabolically within the limits 
of precision, from at least 1 inch on, so that the effects of 
density differences on the flow phenomena must have been 
negligible for this investigation. 

The transverse turbulence levels v' ju behind the grids 
were obtained from the initial rate of spread of the mean 
thermal wake (method due to Schubauer (reference 4)) after 
the effects of molecular spread and finite source bad been 
removed. The resulting levels were somewhat higher than 
those obtained with a hot-wire anemometer but were used 
because of their consistency with the rest of the measured 
diffusion curve. 

Free-stream velocity ffuctuations (without grids) and the 
wake temperature ffuctuations (with grids) were measured 
with the hot-wire anemometry equipment described in refer- 
ence 10. The wires were 0. 00025-inch platinum etched 
from Wollaston; the compensated response of the system 
was flat within ±2 percent over a frequency range from 3 
to 12,000 cycles per second. 

Oscillograms of the temperature fluctuations were recorded 
by photographing a blue oscilloscope tube with fast 35- 
millimeter film in a General Radio type 651-Ab] camera. 

Probability densities of the temperature ffuctuations at 
fixed points in the mean thermal wake were determined 
from photodensitometer traverses of time-exposure photo- 
graphs of a short-persistence (0.001 sec) blue oscilloscope 
tube with the temperature ffuctuations on one pair of plates 
and a 30,000-cycle-per-second sweep on the opposite plates. 
The technique is essentially that used by Simmons and 
Salter (reference 13). 

EXPERIMENTAL RESULTS 

MEAN THERMAL WAKE 

Complete mean-temperature wakes behind a line source 
of heat were measured for 10 different conditions. Arranged 
to indicate the systematic variation of one parameter at a 
time, these conditions were as follows: 


l-in. grid; ^=43.4: U, fps 

8.5 

25.6 

38.0 

C"=25.f) fps; ^=43.4; prid, in 

1 

i 

C’=25.() fps; ’^=86.1; grid, in 

1 

H 

C’=25.6 fps; ,^-in. grid; Jo/Af.- 

43.4 
86. 1 
172.3 

C=25.6 fps; K-in. grid; .t«/A/ 

43.4 
86. 1 
172.3 
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(1)) ^=39, ^=86.1, ^1=^ inch, and 17=25.6 feet per second. 

Figure 3. — Experimental scatter; temperature behind line source of 

heat. 

Here U is tlie mean velocity, M is the grid mesli size, and 
Xo is the heat-source location measured from the grid. Since 
some individual cases enter as elements in two sequences, 
the total number of elements is more than 10. 

Two of the many mean-temperature traverses in the 
^-direction (perpendicular to mean flow and to source line) 
are shown in figure 3 to give an idea of the amount of exper- 
imental scatter. The upper traverse was the worst of the 
lot, even showing an apparent skewness which was not 
borne out by the investigation as a whole. The lower trav- 
erse is more nearly typical of the measured temperature dis- 
tributions from which the standard deviations of the mean 


•8r 


.6 - 


r_ 

M 


A - 


. 2 - 


U 

(ft/sec) 
o 38 
□ 25.6 
O 8.5 


% 

O □ o O 

ft 

O 

o 




_] L_ 

20 30 

Aa 

M 


40 


I 

^ 50 


PTgure 4. — Spread of heat from a line source. 


■^=43.4, M=1 inch. 



Figure 5. — Spread of heat from a line source. ^=43.4, ( ^=25.6 feet 

per second. 


thermal wake were computed. By comparison with the 
reference curve, it is seen to be essentially Gaussian. This 
was the case for temperature profiles at all stations. Since 
the virtually Gaussian character of such a wake has already 
been established by several of the earlier publications, there 
seemed to be no point in reproducing here all of the large 
number of traverses measured. 

The mean-thermal -wake spread for the 10 different con- 
figurations studied is given in figures 4 to 7 as plots of cor- 
rected standard deviation F' against distance from the heat 
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RE f). — Spread of heat from a line source. ^=86.1, f/=25.6 feet 
per second. 



A£ 

M 


Figure 7. — Spread of heat from a line source. ^=172.25, (^=25.6 
feet per second. 


source Ai\ Each point in these figures corresponds to a 
complete transverse temperature traverse. 

In order to have values of transverse turbulence level r'/Z7 
consistent with the thermal-wake behavior, these values 
were determined from the initial angle of spread of the cor- 
rected wake standard deviation (references 2 and 4) instead 
of from direct hot-wire anemometer measurements. The 
results are plotted in figure 8. Since these represent an 
insufficient number of points per grid to permit the drawing 
of reliable curves, some simplifying assumptions were made 
based upon the results of several more-detailed turbulence- 
decay investigations. (See references 14 to 20. The perti- 
nent results of references 14 to 16 are summarized in refer- 


14x10^ 




k2 



(in.) 

xJM 

r 

(fps) 

v'/C 

(porcen t) 

1 

43.4 

8.5 

2.0 

1 

43.4 

25.6 

2.0 

1 

43.4 

38.0 

2.0 


43.4 

25.6 

2.0 

K 

43.4 

25.6 

2.0 


86. 1 

25.6 

1.50 


86. 1 

25.6 

1.40 

y< 

86. 1 

25. 6 

1.30 

V2 

172.3 

25.6 

1.05 


172.3 

25.6 

.95 


Figure 8. — Decay of tur))ulence as determined ))v thermal wake of a 

line source. 

ence 17.) In the light of these papers, it was assumed that 
the decay curves had a common apparent origin and that 
this was obtainable by drawing the best straight line for all 
the available points, independent of the differing wind speeds 
and mesh sizes. In computing any individual Lagrangian 
correlation function from the corresponding wake history, 
the turbulence decay rate was assumed to be given by the 
line drawn through this (‘ommon origin and the specific tur- 
bulence value giving the measured initial spread angle for 
this wake. This is, of course, a very rough procedure, but 
the experimental scatter in tliis whole method of determin- 
ing Eagrangian correlation functions is so great that a more 
extensive study of decay (including the resolution of incon- 
sistencies between wake method and hot-wire method) 
seemed unwarranted at this time. 

TEMPERATURE FLUCTUATIONS 

Distributions of temperature-fluctuation level in the 
thermal wake have been measured by using the hot-wire 
anemometer as a resistance thermometer, that is, at a 
current low enough to render the sensitivity to velocity 
fluctuations negligible compared with the sensitivity to tem- 
perature fluctuations (reference 21). A representative dis- 
tribution of t?'/0 in the a:-direction is given in figure 9. Typ- 
ical transverse distributions are given in figures 10 (a) and 
10 (b). It is clear that the temperature-fluctuation inten- 
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Fir.uRE 9. — Tejiipe rat lire fluctuations at = inches, M=l inch, and 
r = 25.6 feet per second along wake axis. 



(b) Ax =70 inches. 
Fioure 10. — Concluded. 



(a) Ax =10 inches. 

Figure 10. — Temperature fluctuations behind a line heat source at 
-^=43, M=l inch, and U=25.6 feet per second. 

AJ ’ 



of temperature fluctuations. Film speed, 1 foot per second. 

values of Ax with y = 0 . The one-sided and pulse character 
of the instantaneous temperature at a fixed point in space is 
also demonstrated by its probability densit}’. 


sity changes very little with increasing values of Ax. A 
rough explanation of the very high values of (com- 

pared with the concomitant turbulence level, for example) 
in terms of the highly intermittent structure of the thermal 
wake has been given in reference 10 and will be discussed 
in more detail later in this report. This interniittency is 
shown very clearly in figure 11, a series of temperature 
oscillograms recorded at two different positions across the 
thermal wake for a fixed value of Ax and at two difierent 


THEORETICAL CONSIDERATIONS 

HOMOGENEOUS STEADY TURBULENCE AT REST 

For a nondecaying incompressible turbulence with no 
mean motion, Taylor (reference 1) was followed in getting 
an expression for the mean time i-ate of diffusion in the ?/- 
direction (say) from a fixed source as measured by the 
second moment of the probability density of the diffusion, 
that is, the mean-square particle displacement 
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(It dl 


= 2v(l) r v{ti)dli 

Jo 

The bar denotes ensemble average. 

Taking v(l) inside tlie integral, interchanging the processes 
of integration and averaging, and introducing the Lagrangian 
(auto) correlation coefficient (following the fiiiid particle), 
where — /i, 

there results 

y^=2p f rR,(r)dr(/T (1) 

Jo Jo 

This is Taylor’s form. Integration by parts yields a form 
like that in the work of Kampe de Feriet (reference 5): 


correlation function, which he called the ''microscale.’’ 
Applying the same geometrical concept to the present 
function, the Lagrangian microscale 


X 


2 




(5) 


is simply the r-intercept of the vertex-osculating parabola 
of the even function The kinematic significance is 

clearly shown by a series expansion of r(/ + r) in /i\(t): 


K(t) = 


V(t)v(l+r) 




y^=2v^ f‘(i-T)/Mr)dT (2) 

Jo 

In this Lagrangian analysis, v(t) is the velocity of a fluid 
particle in the ?/-direction at time t] v(t — T) is the velocity of 
the same particle at time t — r. Corresponding expressions 
can be written for the rate of diffusion in any direction. 

Diffusion from an infinite line source, the case to be 
discussed here, is a two-dimensional problem in the mean, 
and, in addition to ecjuation (2), 

:Y^=2^ r(t-r)Ru{r)dr (3) 

Jo 

A tensorial generalization of these concepts has been given 
by B atchelor (reference 3) dealing with the behavior of 
XiXj(t) where Xt and Xj are any two of the orthogonal 
displacements of the particle at time t. 

It should be noted that this analysis gives no information 
on the shape of the probability density of F(t) or of X(t). 
In fact there still exists no theory for these. However, 
experiments in flowing turbulence (the case to be considered 
next) show Gaussian density, within the experimental 
precision, for Y{t) at all values of t. 

For this stationary random process, Taylor (reference 1) 
introduced the concept of the Lagrangian "scale,” 

Rt{r)dT (4) 

These have the dimensions of time and are characteristic 
constants of the system. 

In his later work on the (Eulerian) dynamics (reference 2), 
Taylor had occasion to introduce another measure of the 


But Constant, so that 




( 6 ) 


From equations (6) and (5), the Lagrangian microscale for 
v{t) is 


X/ 






or 




X/ 


(7) 


(7a) 


TURBULENCE IN A FLOWING MEDIUM 

The dictates of both practical interest and experimental 
feasibility reqiure anaHsis of the diffusion when there is a 
mean velocity U relative to the source. Since the diffusion 
phenomenon is linear, the probability density (mean-concen- 
tration distribution of tagged particles in the wake) is simply 
proportional to the superimposed probability densities of a 
continuous line of sources moving with the mean velocitv Ij 
with their time (and space) origin at the actual fixed source. 

This is illustrated in figure 12(a) for -^-<1. The circles 

(corresponding to isotropy) are the standard deviations of 
the dispersions that would occur from moving sources. The 
envelope of these circles gives a measure of the mean wake. 
It is obvious that in general the functional form of the mean- 
concentration distribution along a line Ai:=Constant will not 
be the same as the functional form of the same quantity for 
the individual source at time t, that is, at position x=Ut. 
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U 



parabolic behavior, (reference 1), of the diffusive 

T 

process. It will occur for all values of when • "There- 
fore, in this case a simple approximate space-time transforma- 
tion in the mean is permissible, and the variation in Taylor's 
theory of diffusion by continuous movements becomes a 
variation of xlU. (It must be emphasized that the fore- 
going discussion does not apply directly to the possibility of 
applying a space-time transformation to the instantaneous 
turbulence variables. This latter question will be discussed 
later.) 

In the present measurements therefore, the A^-- 

variation of diffusion gives an approximate measure of the 
Lagrangian coz-relation coefficient (in time). Equation (2) 
can be rewritten as: 



(a) ^<1. -V'(0 = -V' F'(0=r' 

(b) 

u 


Figure 12. — Propagation of turbulence from source moving with mean 

velocity U. 


However, it is extremely likely that if (fig. 12(b)) 

the mean concentration along a line Aj-= C onstant becomes 
very nearly of the same functional form as that for the source 
at x=lJt. For a zero-correlation Gaussian density, the 
equivalence is easily demonstrable.^ 


dY 


Tlie condition "^^^1 ''"iT always occur at large enough 

values of A.r when=<[l. This follows from the asymptotic 

U 


3 The most general mathematical restrictions under which the suixii position of a line of 
identical densities will yield a “cross-section” density of the same form have not been studied 
here. It is obvious that the condition of statistical independence is sufficient. The fact that 
the density of each of the velocity components in isotropic turbulence has been found to be 
Gaussian within the experimental precision seems to show that the equivalence under 
discussion is at least a good approximation. 


DECAVING ISOTROPIC TURBULENCE 

When the turbulence is decaying in time (similar to space 
in the flowing turbulence) is no longer constant, and the 
analysis cannot be carried out as far as equation (1). The 
same approach stops with 

( 9 ) 

where the prime denotes root-mean-square value, and 

T> _ V{t)v{t—r) 

^^^-^~v\t)v\t-r) 

There is no a priori reason to believe that fit-r is ^ function 
of r alone, as in the nondecaying turbulent flow. 

At this point Taylor (reference 2) invokes the empirical 
fact that over a wide range of mean velocities (all of which 
give essentially the same distribution of fU in x behind a 
grid) the thermal wake behind a line heat source at fixed x 
appears to be unchanged in form, within the experimental 
error. This is consistent with dependence of the diffusive 
process upon a variable of type 



v'{xi)dxi ( 10 ) 

UJlo 

Therefore, Taylor has postulated the unique dependence of 
Mt-T on the variable rj. With this postulate and the space- 
time transformation valid for small turbulence level, ho 
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arrives at 

(. 1 ) 

or 

( 12 ) 

where 

7 

(lv== (lx 
= v' (It 


Equations (11) and (12) look like the equation for non- 
decaying turbulence. They also give 





(hi 


(13) 


A ph ysical significance of the length rj is underscored by 
the limiting form of ecpiation (13) as rj-^0 and /i\— >1 . Then, 


} =7] 


so that r; is a measure of the lateral diffusion that would 
occur if the lateral velocity fluctuation following a particle 
v(t) remained perfectly correlated but decreased in magni- 
tude according to the decay rate of the turbulence level. 

In postulating Taylor was apparently com- 

paring only diffusive processes in turbulence fields with 

... 

identical (x). Of more general interest is the comparison 

of diffusion in fields with differing turbulence-level distribu- 
tions. Although such a generalized application of his postu- 
late is doubtless not too well applicable, it is conceivable it 
might have approximate success in the more genei’al com- 
parison. For both convenience and lack of any obviously 
superior alternative, his suggestion is therefore applied in 
computing the results of the measurements reported here. ^ 

With the 77-postulate, the /»\ Lagrangian correlation func- 
tion can be obtained from measurements of as a function 
of Ax: 




1 

9 


C; 



(14) 


* After this work was completed, Dr. Batchelor suggested an alternative appro.xiinate 
ai)proach: In order to construct a stationary random function out of the nonstationary v{t) 
one first normalizes the deixmdent variable with its root-mean-square value. (This has been 
automatically accomplished by use of the correlation coefficient.) If the decaying quantity 
is assumed to maintain complete similarity during decay, all characteristic times (e. g., 
Lagrangian time scale and microscale) vary in the same way with t and the new independ- 
ent variable is constructed by dividing t by this f-variation. Unfortunately, this variation 
is unknown a priori, so it would be necessary to assume further that the Lagrangian scales are 
directly proportional to the Eulerian scales, about which there is previous e.xperimental infor- 
mation. II is work has now been published; see reference 22. 


A scale and a microscale can also be defined for Ii%: 



(15) 

(16) 


With nondecaying turbulence, no one-parameter true 
(time) Lagrangian correlation function exists, and the 
77-formulation is much more convenient. A further signifi- 
cance of this variable will appear in the comparison of 
Eulerian and Lagrangian treatments of diffusion from a line 
source in flowing turbulence. 

ACCELERATIONS IN DECAYING TURBULENCE 

A series expansion of v(t-\-T) for decaying turbulence will 
show something about the initial behavior of the true (time) 
Lagrangian correlation function and will indicate an experi- 
mental method for examining a hypothesis of Taylor on the 
interchangeability of instantaneous time and space deriva- 
tives when the turbulence level is low (reference 23). 

Write the Lagrangian correlation coefficient 


Substitute 




v{t)v(t— t) 





(17) 


into numerator and denominator, and restrict the analvsis 
to small values of r: 



Divide numerator and denominator by expand the square 
root in the numerator, and keep terms in 


/fp(^,r) ~ 1 


-[-m. 


^W] 




(19) 


For negligible decay rate this reduces to equation ((>). 

E(iuation (19) shows that a Lagrangian microscale defined 


hy 


1 

W{t) 


liin 

t-K) 


j- l-7f,(<,r) j 


is expressible as 


__rd^E,(h r)l 

" I dp 
x.^(0“2pLv/f/ 


( 20 ) 

( 21 ) 
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III the absence of infoi’inatioii on the algebraic sign of the 
triple correlation term, it is sufficient to require the two 
conditions/ 

1/2 


Introduction of Taylor’s ? 7 -postulate transforms equation 
(19) to 


after the additional approximation that «»'(0) wlien 

7) is very small. 

This gives a new ex])ression for Lagrangian microscale X,: 


J__J_ r t— Yl 

V~2’pLv/i?/ \dr?/_l 


(23) 


Equation (23) is in contradiction to equation (17) of part 
IV of reference 2. In that work Taylor has apparently 
assumed that v' is a nondecaying function of i\. However it 
certainly is decaying, even in teims of this distorted coordi- 
nate, and the dc'/dij-term must he included. 

Since X, and v'{x) can he detennined experimentally from 
the mean thermal wake behind a line heat source, equation 

(23) permits determination of (^) "'hich is simpl}’ related 

to the mean-square “Stokes” acceleration (^^) • 

This quantity is of particular interest for the possiliility of 
an instantaneous space-time transformation at low turbu- 
lence levels. This was first proposed 1)3" Ta\"lor (leference 23) 
and has since been used very widely, especially to get approx- 
imate values of partial derivatives with respect to x (the 
mean-flow direction) by measurements of time partial 
derivatives. 

The total (or Stokes) derivative of v{x, y, z, t) in a turbulent 
ilow with mean velocit}^ U along the ^-direction is 


dv dv , .TJ,. 

*-5T+<‘'+“' dl+' 5S+*’ 51 

'I'aylor’s hypothesis amounts to the statement that 


(24) 



dv ^ jjdv 
dt ~ dx 


or 

dv jT dv , 

(24a) 

with 

6^ 

1/2 

<1 

(25) 


—2 (bv\ 

w J 


In detail, equation (25) is 

/'dvV . bv 




dv dv 
dXf dxj dt dxj 




LI 


\bxj 


1/2 


<1 


(25a) 




L-’(i^yj 


UiU) 


dv dv 

dxi dXj 


<1 


1/2 


(26) 


L^’(S)’J 

]5„t (^y fan be determined from measurements of Y'{x); 


( 


= 2%, where X is the Eulerian microscale (reference 2) ; 

,dx/ X^ 

and upper bounds, in terms of measurable functions, can he 

set on UiUj ^ ^ with the use of Schwarz inequalities. Thus, 
oXi ox j 

an experimental check of the requirements in equation (26) 
is to be made in the section entitled ‘^Computation of 
Results. 

RELATION BETWEEN EULERIAN AND LAGRANGIAN MICROSCALES 

Taylor (reference 2) inferred an approximate relation 
between and X by neglecting the effect of viscosity on 

pressure-gradient fluctuations and estimating J as 

being ap])roximatel\" 3p ^ constant 

ratio \J\ for all turbulence. The rough nature of this 
anal 3 ^sis induced Heisenberg (reference 25) to conduct a 
more detailed study of the static-pressure fluctuations and 
to reestiniate the \J\ ratio. However, he followed Taylor 
in ignoring the dv' /drj-Un'm in the relation between X,, and 

equation (23)) and in neglecting viscous terms in 

/ dp\^ 

the relation between \ and * 

Although these omissions are probabl}" not serious except 
in the low Re 3 'nolds number range, it seems interesting, if 
onl}^ for the sake of completeness, to use a Heisenberg type 

of approximation for repeat his treatment with 

the omissions rectified. 


5 Lin has discussed the validity of Taylor’s hypothesis using a slightly different formulation 

in reference 24. He points out there that if = L <^hen (from the Schwarz inequality 
^ i 

2 ''■here a, is a set of numbers. 
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l^rom the complete Navier-Stokes equations in the wave- 
number space, Heisenberg deduced an approximate expres- 
sion for (Vp)2 in terms of mean quadruple products of the 
harmonics'' of the velocity field. His principal simplifying 
assumptions were; 

(a) Different Fourier components of the velocity held are 
uncorrelated 

(b) The turbulent energy spectrum is given by the solution 
to his equilibrium-energy-transfer equation, above a lower 
cut-off wave number ko 

Following these, but using Chandrasekhar's (reference 26) 
solution to the Heisenberg equation instead of the inter- 
polation formula used by Heisenberg, there results 


( 


^2// 


25.2 


P ,^3 


'^3 


(27) 


where k is a dimensionless empirical constant. 

The numerical constant in equation (27) would perhaps 
have been given more accurately by the use of a ''self- 
preserving" spectrum (calculated by Chandrasekhar from 
Heisenberg's equation) instead of the stationary spectrum 
with low cut-off wave number. Time was not taken to make 
the requisite additional calculations because: (a) The value 

principally upon the high-wave-number 

region of the velocity spectrum rather than the low-wave- 
number region, where the difference would be greatest, and 
(b) the experimental results and (especially) the value of 
K both have a considerable range of uncertaint 3 L 

The mean square of the y-component of the Navier-Stokes 
equation will lead to a relation between and X: 

dv dv . dv , dv , dv 
(It bt dx dy dz 


therefon 


1 > T7'> 

= ^+pV‘-v 

P 




(28) 


where the correlation between pressure gradient and v(‘locit 3 " 
J^aplacian function is zero because of isotropy. 

Equations (23) and (27) give the first two terms in terms 
of the microscales, and the mean-square Laplacian function 
is expressible in terms of the fourth derivative of the Von 
Karman-Howarth /(r) correlation coefficient at r=() (refer- 
ence 27): 


Qr: 

(V^vr=^ {vrp\o) 


(29) 


Consequentlv, equation (28). becomes 

The second term in tliis equation can be replaced by the 


t urbulence-decay equation 


d{vy_ 

dt 

whence 

dv' 
dt ~ 

so that 



{v'f 


(30) 




W(0) 


therefore, 




K « ft ft’ ' 6 ■' ""(?? 


2 \ 2 


(31) 


wliere R\ 


v'\ 

V 


Batchelor and Townsend (reference 19) have deduced an 
expression for /^*’(0) which is valid in the region of decay 
wliere both l/(y')^ and X^ increase linearly with t (correspond- 
ing to large values of 7?x): 


XV‘'(0)=^+|t?xS 


(32) 


ewness factor. Their 


where S=-(|j)Y[(gy] , ,1„. sk 

experimental results showed tS=0.39 approximately constant 
for isotropic turbulence. Then the estimate for X/X, 
becomes 




(33) 


I he value of k was first estimated by Heisenberg (reference 
25), from measurements of turbulence decay, as 0.85. 
This method may be regarded as emphasizing the (relativelv 
low wave number) energ^^-bearing range of the spectrum. 
Lee (reference 28) worked out an estimate based upon 
skewness factor (/c=0.13), which gives heavv weight to the 
high-wave-number range. Proudman (reference 29) has 
reestimated k bv comparison with measured curves of the 
double and triple velocity correlations. The value k=0.45 
leads to reasonably good agreement for the moderately 
high-wave-number region, over a wide range of values of R\. 

It may be remarked that the supposed constancY" of k is 
merelY" a postulate of the Heisenberg dimensional formulation 
of the spectral transfer function. In fact it is by no means 
obvious that this turbulent part of the transfer is quantita- 
tivel}' independent of the amount of s])ectralh^ local dis- 
sipation to heat. In anv case, Proudman's estimate of 
/c=0.45 has been used here. Therefore, 


(0 


12.5 29 


(34) 


In the limit of R\ 0, equation (34) does not apply since 
equation (32) does not apply. However, the appropriate 
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limiting relation ean be obtained directly. In this limiting 
condition the pressure term in equation (28) is negligible 
compared with the viscous term (the former goes with 
l//?x and the latter with l//?x^), and the Eulerian velocity 
correlation coefficient is (reference 27) This 

gives /'(O) =3/X^ and 



(^-direction). Of course, &a has been neglected relative to 
QU in equation (39). 

Equation (3*8) has two unknowns, and the first objective 
is to express^ as a function of the (more easily measurable) 
'0{x,y). After integration with respect to y, 

'^=-uj“^dy + F{x) 


EULERIAN ANALYSIS OF HEAT DIFFUSION FROM A LINE SOURCE 

'Phe two-dimensional turbulent-heat-transfer equation is 


dx oy pCp \02 


Q_^d 
dx'^ dy 




(36) 


where 0 is mean tempei'ature, is temperature fluctuation 
about the mean, k is thermal conductivity, and Cp is specific 
heat at constant pressure. 

For the thermal wake behind a line source_in isotropic 
turbulence with constant mean velocity F=0. With 
restriction to low turbulence level, a ^0)oundarv-layer^^ type 
of approximation can be applied to the mean wake, so that 


and 


ye 

dy^ 


d /— N 


so that oquatioii (36) takes the a])])roximate form 



k ^ 


C)7J 


(t^v) 


(37) 


But, by symmetry, t?i»=0 for y=0. Hence F(x) = () and 

This relation is sufficient for the computation of &v{x,y) 
from the measured B(x,y)_hut the empirical fact of simple 
geometrical similarity in 0(x,t/) suggests exploitation of the 
consequent simplification. 

Assume 

e(x,y) = eo{x)f(^) (41) 

where ^=i//F'(x). This transforms equation (39) to 

OXx)Y'(x)=-^ (42) 

where H*= =Constant. It transforms equa- 

2pc,Jo /(Qrf^ 

tion (40) to 

(43) 


It must be emphasized that for this particular initial con- 
dition on the temperature (effectively a “]^int source”), 
the restriction to small turbulence level v'/U<^l does jiot 
imply that is small. In fact, for this problem is 
often greater than unity, es])ecially at the '‘edge” of the mean 
wake, as has been discussed in reference 10 and will be 
brought out again later in the present report. 

When the molecular transport can be neglected relative to 
turbulent transport 


With t'quation (42), 0^(x) can be eliminated from equation 
(43), and after integration of the dimensionless integrals this 
leads to the final form for the turbulent-heat-transfer 
correlation, 

M., y)=H-UW ( 44 ) 

where 

^ Y'(x) 


j-7 50 5 

5x by 


(38) 


an equation given in reference 10; a slightly more general 
treatment follows: 

With a constant rate of heat generation (similar to steady 
state in the average), the application of a Von Karman 
integral-relation treatment to equation (38) yields an integral 
condition: 

^ CO 

2pCpU \ 0 rf// = Constant=// (39) 


where // is the average time rate at which heat^crosses all 
planes perpendicular to U per unit length of heat source 


The same sort of analysis can be made on equation (37) 
which includes the molecular conduction, but the rather large 
experimental scatter in the present measurements seems to 
make such a refinement inappropriate. 

An "exchange” coefficient or "diffusion” coefficient for 
turbulent heat transfer kr is simply expressible in terms of 
0(x, y). A conventional procedure for semiempirical analyses 
is to write for the turbulent transport an expression just like 
that for the molecular transport: 


— pCp^v = kr 


by 


which serves as the definition of 


( 45 ) 
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For the simple case of equation (44), it turns out that 


kj' — pCpU 


■f/(0 

df 

. J 


(ly 


(46) 


This has the particularly interesting property that in a 
nondecaying turbulence at very large values of where the 
mean thermal wake spreads parabolically (reference 1), kr 
becomes independent of explicit dependence on x. 

A more startling simplification follows for a particular 
mean temperature distribution across the wake: All depend- 
ence of kr on y disappears if 


that is, if 




.m = e 


-hv- 


(47) 


But this is the Gaussian function, which is found emi)irically 
to fit all the measurements within the experimental scatter. 
Hence one arrives at the empirical result that in both non- 
decaying and decaying turbulence kr is independent of y 
in the thermal wake behind a line source of heat. From 

equation (47), 


kr= pCpUY'{x) 




dx 


(48) 


It can be seen that in the nondecaying case at very large 
values of X, kr is constant and independent of both y and x. 

RELATION BETWEEN SOME LAGRANGIAN AND EULERIAN 
PARAMETERS IN TRANSPORT 

There has apparently been little effort to relate the 
Eulerian and Lagrangian formulations of turbident diffusion 
up to the present time. Exceedingly simple boundary con- 
ditions permit some connection to be made in restricted 
ranges of the present problem. 

For nondecaying or (with inTich less accurac}^ decaying 
turbulence, equation (13) ap])Iies: 


Y^ = 


= {r] — rji)Rr,{vi)dvi 


or, equation (12), 


dr] dr] 


= 2 Jo li^{r]i)dvi 


Tj ^ 

For small values of 77, Ryj= 1 — whence 


('-is?) 


(49) 


drt \ 3 X,v 

Substituted into equation (44) these give 

while equation (46) becomes 

2. m /j 1 

wliere 2 ,=l Jo%'* ami In fact, 

V ^ 

for small values of rj, n «= x. 

U 

v' 

For losv-level nomleeaying turbulence, t\==x, and equa- 
tion (51) becomes 


(50) 


(51) 


(52) 


§v=iri{o^J 1 -— ) -- 
12r’x,V2;' 

while equation (52) becomes 

5 ill 

<(/<(£( 4(-X?/ 

As 77->0 l)oth equations (51) and (53) reduce to 

while equations (52) and (54) reduce to 
kr=pc,v' y 

At the other cxti'eme, svlien y is very large, 

T^=2„L,-2.V, 

n 00 

where 4/i = J 72/?,(t;) d») = Constant and 

(ly 

In this case. 


(53) 


(54) 


(55) 


(56) 


(57) 


, =2i, 

(In 


dv=Il*m) 


yL/= 


and 


[2(nL,-.V/,)l^^^ 


kr = pCpV' 




(58) 


(59) 


(60) 


djjdi [2(2,L,-A/,)]''^ 

As indicated previously, the Tj-variation can be expressed 
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in terms of x or 

If j) is allowed to become large enough to make 


v' 


— yy 

(61) 


(62) 


The above fonnulas take on particularly simple forms if 
the empirical result of a Gaussian /(^) is utilized: 




Then the general expressions for dr and kr (equations (44) 
and (46)) become 


dv{x,y) 


H* dY' 
{Y'f{x) dx 



(63) 


and 

f]Y' 

kr(x,y) = pc,Y'{x)^U (64) 


the latter having been deduced in the previous section. The 
])articular forms for small values of ij would follow from sub- 
stitution of equations (49) and (50) into these two. 

However, the most interesting form occurs for very large 
values of v- There Y'ix) is given by equation (57) and 



_ n*yL,^ 

d»(x, ?/)~ [2 (t,Z,-Mj)]5/' 

y) = pCpV' L^=Const2i\\t 
For still larger values of rj, such that Mi<^t}L, 


(65) 

( 66 ) 


and 


H* 

^v{x,y)=-^ 


U 


2v'2 


1/2 



kf[x,y) = pCpv'Lrj = Constant 


(67) 


The constanc}^ of kr for large values of rj (large values of t 
or x) is to be expected; a treatment of molecular diffusion by 
this method must certainly yield a constant coefficient for 
times much larger than the mean free time (of flight) of the 
molecules — that is, for all “macroscopic times/’ Put an- 
other way, the simple parabolic beha^dor of Y' for large 
values of r; is a sure indication that 0 obeys the simple 
classical diffusion equation with constant coeflicient, when 
viewed extremely “coarsely.” 

Perhaps the chief interest of equation (66) is its identifica- 
tion of as a significant Lagrangian length for diffusion at 
a large distance from the source. It enters the expression 


for turbulent-diffusion coefficient in much the same way as 
mean free path enters the expressions for the molecular- 
diffusion coefficients. Furthermore, its role appears to be 
much like that attributed to PrandtPs “mixing length,” 
which was brought into the turbulent-transport problem in 
a more or less intuitive fashion. 

Of course, the possible crude nature of Taylor’s original 
r;-postulate may render the significance of Zr,more qualitative 
than quantitative in the case of decaying turbulence. 

COMPUTATION OF RESULTS 

Although Taylor's assumption of the unique dependence 
of upon 7 ] is not likely to be accurate for collapsing 

together cases with widely differing turbulence decay rates, 
it does provide a relatively simple relation between F'(x) 
and Therefore all of the mean-thermal-wake data 

were reduced on the 17-basis. 

In principle the complete A\(i7) curve can be obtained from 
T^(x) by double differentiation (equation (12)): 


or 


R (F"^) 

R =i(^— ^ — (F^) 


( 68 ) 


However, simple double differentiation of the squares of 
a curve as uncertain as F'(x) seems almost hopelessly in- 
determinate— although Taylor (reference 2) and Collis 
(reference 7) have apparently followed this procedure. A 
somewhat more circumspect technique has been tried here: 
The values of X, and Ly, were determined first, through certain 
limit relations (to be described). Then the curve was 
determined b\' double differentiation, subject to the restric- 
tions of agreement with the previously determined scales. 
Thanks to rather poor determinancy of values of and 
this method is not so much of an improvement as it might 
first appear. 

LAGRANGIAN MICROSCALE X, 

If equation (13) is restricted to very small values of 17 the 
parabolic approximation for can be introduced: 


therefore 


whence 


F^=2j\,-,n)(l-^) dv^ 



(69) 


(70) 


The computational procedure was to plot FV^;^ against 17^ 
and to estimate the slope of the faired curve at rj = 0 where 
the curve must pass through unity. The abscissa intercept 
of the 0-tangent is 6X/. The actual points, faired curves, 
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(a) X, = 0.36. 

(b) L, = 0.43. 

Figure 13. — Determination of X, and L,;^=43.4, M=\ inch, and 
U=8.5 feet per second. 


and tangents for all cases are presented in figures 13 to 22. 
Clearly the precision is jioor. 

LAGRANGIAN SCALE L, 

Consider equation (12) in the limit as rj-^oo. It immedi- 
ately gives 



and the graphical procedure based on this is also presented in 
figures 13 to 22. Some of the asymptotic slopes drawn are 
not the best representation of the experimental points. This 
is due to the auxiliary (assumed) restriction that /?, cannot 
increase with increasing values of rj as long as has not 
previously" dropped below zero. The graphical jirecision 
attainable is perhaps a little better here than that for X/, but 




(a) X,=0.34. 

(b) L, = 0.33. 

Figure 14.— Determination of X, and Lp; ^=43.4, A/=l inch, and 
(7=25.6 feet per second. 


the square root necessary to get X, means that X, is determined 
about as well as is i,. 

LAGRANGIAN CORRELATION FUNCTION 

AVith X, and determined, the initial (small Ax) and final 
(large Ax) behavior of the curve F'(Ax) is prescribed. These 
parts of the curve were drawn on a graph with the experi- 
mental points. Then the fairing in of a reasonable cen- 
tral portion to this mean F'(Ax) curve was a relatively 
simple matter. The curve was then obtained by double 
differentiation. 

The curves drawn for Y' {Ax) in figures 23 to 32 were de- 
termined in the fashion described above, as were the curves 
for 7?^ in the same figures. 

EULERIAN MICROSCALE X 

In view of the approximate nature of the determination of 
no new direct measurements were made of the Euler ian 
microscale X. Instead, X w-as computed with the energv 
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(a) X, = 0.27. 

(b) L, = 0.33. 

Figure 15. — Determination of X, and L,; ^=43.4, 3/=! inch, and 
r^=38 feet per second. 




(a) X,=:0.25. 

(b) L,=0.34. 

X 1 

Figure Ih. — Determination of X, and L,; ~=43.4, and 

(7=25.6 feet per second. 


equation for isotropic turbiiloiice, from the measurements of 
turbulence decay: 


X2=-10r-J:, - 
dv^ldt 


(72) 


or, with the space-time transformation, 


X2=_i0i. — 

uT 

dx 


(73) 


at the California Institute of Technology (reference 18) on 
grids of essentially the same geometry. 

Table I summarizes the results for Lagrangian and Eulerian 
scales and microscales. The results have been grouped to 
show the effect of systematic variation of one parameter at a 
time. Some of the results are presented in figures 33, 34, 
and 35. 


INSTANTANEOUS SPACE- TIME TRANSFORMATION 


EULERIAN SCALE L 

Earlier investigations have shown that the Eulerian scale 
in a grid-produced turbulence is closely a linear function of 
the mesh size of the grid producing the turbulence (for a 
given value of x and grid geometr}^ and is not significantly 
dependent upon the mean velocity (or grid Reynolds num- 
ber, provided it is sufficient to cause turbulence) . Therefore 
the values of L have been deduced from earlier measurements 


The permissibility of an instantaneous space-time trans 
formation in flowing turbulence. 


^ ^ 

^x^ U dt 


(74) 


can be estimated in accordance with equations (26). For 
equation (74) to be valid, the sufficient requirements are those 
given in equations (26) that: 
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(a) X,=0.24. 

(b) L, = 0.25. 

X 1 

Figure 19. — Determination of X, and L,; 

I7=25.() feet per second. 

value. For isotropic turbulence; equation (76) can be 
written 

dv dr ^ ^ r^\ {v^ 




(77) 


'riuis, the second condition in equation (26) will be satis- 
fied if 


2.7 =<C1 

U 


(78) 


Both T and r'/T7 for tlie flows studied are presented in 
table I. It is clear that for these flows instantaneous x and t 
partial derivatives may be taken proportional with reasonable 
confidence. 




(a) X,=0.1I. 

(b) L,=0.19. 

Figure 20. — Determination of \ and L,; b inch, and 

t/=25.0 feet per second. 


f?i?-C0RRELAT10N 

ddie Eulerian measure of transverse turbulent heat trans- 
port is computed from the mean temperature distribution, 
ddie dimensionless form, &vlOoU, is given for two typical cross 
sections in figures 36 and 37. 

The measurements of v'/U and of r?'/0o permit calcula- 
tion of the correlation coefficient and this is also 

given in figures 36 and 37. 

For the convenient and reasonably accurate assumption 
of Gaussian mean temperature distribution, the correspond- 
ing turbulent-heat-transfer coefficient kr follows from equa- 
tion (64). It was found to be independent of ?/, and typical 
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(a) X, = 0.22. 

(b) L, = 0.22. 

1‘Tovre 21.— Delorminatioii of X, and L,; i^=l72.3, inch, and 

U = 2o.i\ feet per second. 

curves of /tt/Ic are given in figure 38. 'J'he data for kr at 
three difl’erent speeds behind the 1-ineli grid are roughly 
collapsed together through division of kr by pCpV'L,,, as 
suggested by ecpiation (66a), an asymptotic result foi- 
nondecaying turbulence (ng. 39). 

DISCUSSION 
LAGRANGIAN VARIABLES 

h>ven a cursory examination of the technique used in this 
investigation for the determination of Lagrangian correla- 
tion shows that, as physieal measurements go, this method 
is a “had” one, largeh’ because of the inherent double 




(a) X, = 0.17. 

(b) L,=0.13. 

Figure 22.— Determination of X, and L,; ^=172.3, -V=|- inch, and 
f/=25.6 feet per second. 

differentiation between measured variable and desired 
information. 

higures 13 to 22 suggest an uncertainty in values of X, 
and Z,, as large as ±20 percent, in spite of moderately good 
precision in the measurement of individual temperature 
distrihutions such as the lower curve in figure 3. 

As mentioned earlier, the values of v' fU eomputed from 
initial wake spread are consistently higher than those 
measured with the hot-wire anemometer (reference 18). 
The same relative result was encountered during a brief 
investigation following that reported in reference 18. Up 
to the present time there has been no satisfaetory explana- 
tion of the discrepancy. A tentative hypothesis which would 
at least aceount for its direction may be based upon a human 
weakness in the visual averaging of the reading of a fluetua- 
I tion pointer; there seems to be a tendency to choose an 
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(a) Spread of heat from a line source. 

(b) Correlation function 

Figure 23. — Spread of heat from a line source and correlation function 
Rrj for M=1 inch, f/=8.5 feet per second, and ^=43.4. 

^Siverage’^ more or less halfway between the extremes of the 
needle travel. Thus a pointer motion with very skew prob- 
ability density (greater than 0) would tend to be averaged’^ 
at too high a value. The thermocouple voltage in one of 
these thermal wake traverses has just this character (fig. 11). 
Hence a visual averaging might yield too high a wake width. 
If this effect is nonnegligible, it is advisable to employ some 
electrical means of averaging for skew signals, for example, 
the ffuxmeter and bucking circuit described in reference 10. 

In view of tlie considei-able uncertaint}^ in as well as 
that in X and v' the poor degree of agreement between 
experiment and theory shown in figure 33 is understandable. 
Since the two undetermined constants in the theoretical 
result have been evaluated from sets of o'xperiments com- 
pletely independent of the present ones, this agreement can 
be viewed as an affirmative result. 

Since some sort of Lagrangian scale should be a significant 




(a) Spread of heat from a line source. 

(b) Correlation function 

Figure 24. — Spread of heat from a line source and correlation function 
for M=l inch, t/=25.0 feet per second, and ^=43.4. 

length in turbulent heat and mass transport, as demonstrated 
in the analytical section of this report, an effort has been 
made to find some systematic variation in the values of 
Figure 34 might be construed to indicate a monotonic 
decrease of LJL with increasing values o( It is inter- 
esting to note that a decrease was also observed for the 
ratio of mixing length to tube radius by Nikuradse (reference 
30) in fully developed turbulent tube flows. In order to 
determine whether these two rates of decrease with increasing 
Reynolds numbers are of the same order of magnitude, an 
estimate has been made of the magnitudes of Rl correspond- 
ing to Nikuradse^s results given in figures 28 and 29 of 
reference 30. Both scale-to-diameter ratio and average 
turbulent levels for various Reynolds numbers were esti- 
mated with the help of Laufer’s data on turbulent channel 
flow (reference 31) at various Reynolds numbers. The 
absolute level (i. e., the ordinate scale) of the resulting ImaxIL 
against Rl curve was adjusted to give the most reasonable- 
looking fit with the LrjIL data. This is the dashed line in 
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(a) Spread of heat from a line source. 

(b) Correlation function ’ 

Figure 25. — Spread of heat from a line source and correlation function 
Rrt for 3/=l inch, Lf=38 feet per second, and -^=43.4. 

figure 34 and it shows at least a qualitative resemblance. 
It is likely, however, that LJL is not a unique function of 
Rl. 

Of course not all of the scatter in figures 34 and 35 (involv- 
ing can be attributed to simple lack of experimental 
precision. Some is evidence of the fact that the Taylor 
postulate of Lagrangian correlation function being uniquely 
a function of rj is certainly not very closely true. Further- 
more, taV)le 1 does show rather systematic variations of L,, 
in some of the three-point groups. Most noticeably, there 
is a regular decrease in Z,, with increasing Xo/M (or perha])s 
with decreasing v'fU) for each of the three grids. 

TEMPERATURE-FLUCTUATION FIELD 

Fairly close behind a line heat source in turbulent flow, 
the random pulse nature of the temperature fluctuations at a 
fixed point has been established ])V the oscillogram in refer- 
ence 10. This is confirmed by the first two oscillograms in 
figure 11, with their highly skew probability densities at 




(a) Spread of heat from a line source. 

(b) Correlation function Ry,. 

Figure 26. — Spread of heat from a line source and correlation function 

1 — T 

Rr, for M=— inch, t/=25.6 feet per second, and -^=43.4. 

A2’/il/=10. One of the objectives of the present investi- 
gation was to find out whether this distinctly pulsed char- 
acter persisted far downstream or whether molecular heat 
conduction becomes increasingly effective in smearing out 
tlie pulses, until they are no longer distinguishable as such. 
The third oscillogram and probability density in figure 11 
(Ax/M=70) does show^ a decided trend away from the pulse- 
type signal. The molecular broadening of the laminar w'ake 
(corresponding to the pulses) decreases the relative si)acing 
of the pulses in i?(/) at any point in the ‘Turbulent wuvke’’ 
region. This is a reduction in relative length of the flat 
(?^=0) base lines between pulses, giving greater statistical 
symmetry in about its mean, that is, reducing the skew^- 
ness of 

A simple analysis w ill show the existence of an asymptotic 
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(a) Spread of heat from a line source. 

(b) Correlation function Rr,. 

Figure 27. — Spread of heat from a line source and correlation function 
Rr, for inch, 17=25.6 feet per second, and ^=43.4. 




(a) Spread of heat from a line source. 

(b) Correlation function Rr,. 

Figure 28. — Spread of heat from a line source and correlation function 
Rr, for A/=l inch, l7=25.6 feet per second, and ^=86.1. 


behavior of molecular-conduction effects in a nomlecaying 
turbulence. For a nondecaying flow turbulence and very 
large values of t{-xlU) the mean-square wake spread due 
to turbulent motion is 

(79) 

On the other hand, an approach to molecular diffusion 
through Taylor’s eonce])t of “continuous movements” gives, 
for any macroscopic distance downstream, the mean-square 
thermal wake width, 

1V«2A^x (80) 

where A is the mean free path and c is the root-mean -square 


molecular velocity. 

From equations (79) and (80) 



For a typical case, take L=l centimeter, ?;' = 10 centi- 
meters per second, A = 6X10“® centimeter, and c=5X10^ 

Y ' 

centimeters per second. Then 0.1 7. For people accus- 

tomed to thinking of molecular transport as negligibly small 
in turbulent flow (e. g., in shear flow), this ratio will appear 
quite large. The values of kr/k plotted in figure 38 also 
show that at these low turbulence levels the molecular 
thermal conductivity is not necessarily negligible comiiared 
with the tuz-bulent transport. 
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(a) Spread of heat from a line source. 

(b) Correlation function /f,. 

Figure 29. — Spread of heat from a line source and correlation function 

1 — X 

Rr, for M=~ inch, U=2o.6 feet per second, and -^=86.1. 

The tomperatiire-fluctuation-lev(4 distrihutions d'/O across 
the wake (figs. 10 (a) and 10 (h)) show the same character 
as tliat measured at much higher turbulence level in a jet 
(reference 10), with somewhat lower minimum values, which 
are attributable to the lower turbulence level. A rough 
evaluation of the behavior of the statistical variables in this 
turbulent thennal wake is obtainable by recalling that it 
consists of a randomly ‘Svaving” laminar thennal wake. If 
t?(0 is crudely represented by a randomly spaced sequence 
of identical rectangular pulses with height width j, and 
average spacing .9, it is easily seen that 





(a) Spread of heat from a line source. 

(b) Correlation function 7?,. 

Figure 30. — Spread of heat from a line source and correlation function 
Rn for M=— inch, f/=25.6 feet per second, and ~=86.1. 

of the turbulent thennal wake have higher values of sjj, the 
behavior of equation (83) is consistent with the experimental 
distribution. If the analysis were repeated with triangular 
pulses, for example, the quantitative estimate would doubt- 
less be more realistic. The higher values of (^'lo)min en- 
countered at higher values of v'jU are an indication that 
for a given width of laminar thermal wake, the higher value 
of v'lU leads to a higher minimum value of sjj. 

The closely Gaussian shape of e/e^ against y has already 
been pointed out. If 



be introduced, there results 


This permits ??Vb to vai v between 0 and 00 as sjj travels the 
permissible range from I to 00 . Since points nearer the edge 



( 84 ) 
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(a) Spread of heat from a line source. 

(b) Correlation function 

Figure 31. — Spread of heat from a line source and correlation function 
Un for inch, U =25.(5 feet per second, and ^=172.3. 


and 



Both of these expressions liave behavior consistent with the 
experiments. 

The form of dimensionless transverse turl)ident-heat- 
transfer rate can be deduced for small values of Ax 

(such that with this pulse representation of In 

this picture, is the correlation between a continuous ran- 
dom variable v and a random pulse signal -d- which “fires’" 
every time the continuous variable passes through a specific 
value 


Therefore 

s Ax 




(a) Spread of heat from a line source. 

(b) Correlation function R^. 

Figure 32. — Spread of heat from a line source and correlation function 
Rr, for inch, U =25.(5 feet per second, and ~= 172.3. 


dv 



( 86 ) 


where s and 0 are functions of y. 


With Gaussian 0(?/), 


oT^ 




(87) 


ddie direct comparison between this crude picture and the 
experimental results will be confined to the correlation co- 
efficient = This is of particular interest in view 

of the surprisingly higli experimental values. With equa- 
tion (85) and the fact that v' fU=Y' jAx, there results 





( 88 ) 


This contains the undetermined constant which can 

be obtained from any one of several experimental results. 
Figure 36 includes one plot of equation (88) with de- 
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Figure 33. — Ratio of Lagrangian and Eulerian microscales. 
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Figure 34. — Ratio of Lagrangian and Eulerian scales. 

termiiiud from (Hiiialion (84) and the experimental value of 
\^'lo)min and one plot with determined by matehing 
e(|uation (88) with the experimental result at y/Y'=lA. 

It is also surprising to find that the experimental (R^v)max 
at large values of Ax is even larger than that at small values 
of Ax. 'Fhis may be due to a eonsiderable experimental 
error; the resistance-thermometer voltage signal is much 
lower here. Ihifortunately no relation corresponding to 
equation (88) has been deduced for large values of Ax, 
where /?, is essentially zero. 

The criteria for Taylor’s hypothesis of the interchange- 

abilitv of instantaneous space U ^ and time derivatives 

dx dt 

(assumed by him to depend only upon turliulence level) 
have been expressed in equations (75) and (78) as functions 
of turbulence level and X/X,,. If X/X^ is replaced by its 
theoretical expression (equation (84)) in terms of H,,, 
equation (75) becomes 
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Figure 35. — Ratio of Lagrangian"microscale and scale. 
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(b) against y/Y'. Dashed curves are result of analysis based on 
rectangular temperature pulses. 

Figure 36. — Heat -transfer correlation across thermal wake, computed 
from measured mean temperature distribution. 11 = 25.6 feet per 
second, M= \ inch, Xo=43.4 inches, and Ax = 10.5 inches. 
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, . i)v t y ■ 

(a) against ^ 

ejj ^ 

(t>) /?,?, against ylY'. 

Figure 37. — Heat-tran.sfer correlation across therinal_ 2 'ake, computed 
from measured mean temiierature distribution. (7=25.6 feet per 
second, A/=l inch, x„ = 43.4 inches, and Ax = 70 inches. 


U 


Rj 


<1 


(89) 


8x10 




(a) Turbulent-heat -transfer coefficient at three airspeeds. T/=l inch 

and |j=43.4. 

(b) Turbulent-heat-transfer coefficient for two different grids. 

I/ = 25.6 feet per second and ^=43.4. 

Figure 38 .— Turbulent-heat-transfer coefficient for three airspeeds 
and two different grids. 


For estimates of most flows tlie first term in the parentheses 
ean be neglected; values of R\ less than 5 or 10 are rare. 

vSince equation (34) has now been rouglily verified by 
experiment, equations (89) and (78) may serve as approxi- 
mate criteria for the validity of Taylor^s hypothesis. 

In the limit of R\-^0 when equation (35) replaces equation 
(34), there follows a simpler criterion to replace equation (75) : 

±(fy< m 

SUMMARY OF RESULTS 

The following results were obtained from the investiga- 
tion of the difltusion of heat from a line source in isotropic 
turbulence. 



Figure 39. — Dimensionless turbulent-heat-transfer coefficient. 


■^=43.4 and A/=l inch. 
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1. I he thermal wake behind a heated wire set perpendie- 
ular to a flowing isotropic turbulence (at sufficiently low 
wire Reynolds number) consists of a randomly “waving” 
thin, laminar, thermal wake wliose variations in lateral posi- 
tion give what may be called the turbulent thermal wake. 
At a fixed point not too far behind the wire the instantane- 
ous temperature difference is a random pulse function, 
and the nature of the turbulent heat transfei- can be deduced 
on this basis. Farther downstream the distinct pulse nature 
tends to disajjpear. 

2. The mean transverse temperature distribution £- (y) 

appears to be Gaussian within the experimental jjrecision 
for all distances behind the wire. 

3. An Eulerian analysis of this turbulent-heat-transfer 
inoblem permits computation of the turbulent-beat-transfer 
eoefficient kr which is essentially constant with respect to 
the distance in the direction of the measured diffusion ?/ for 
these boundaiy conditions. It is found that at low turbu- 
lence levels (approximately equal to 1 to 2 percent) the 
molecular beat transport is not vanishingly small compared 
with the turbulent heat transport. 

4. Although Taylor’s postulate that Lagrangian correla- 
tions in deca\-ing turbulence can be made similar by intro- 
duction of an independent variable i: v' (t)dl (wliere, t is 

time and v' is the root-mcan-square instantaneous veloc- 
ity fluctuation in the y-direction) seems to be an oversim- 
plification, it has been applied here for convenience in the 
reduction of data. A simple comparison of Eulerian and 
Lagrangian analyses for diffusion in nondecaying turbulence 
shows that for large values of the distance from the beat 
source Aa: the Lagrangian scale Z, entei-s the expression for 
kr, the turbulent-heat-transfer coefficient, much like the 
empirical mixing length in the old turbulent transport the- 
ories. Therefore some properly modified generalization of 
Ta 3 ’lor’s jj-postulate should prove useful. 

5. A correction and generalization of Heisenberg’s theo- 
retical expression for the ratio of Eulerian to Lagrangian 
microscale X/X, as a function only of the turbulence Rey- 
nolds number based on microscale /?x has been made and 
seems to agree roughly- with experiment. It must be noted 
that since X, depends only upon a transformation = 
and not upon the integral postulate stated above, its valid- 
ity is not impaired by any failure of the. integral postulate 

6. Taylor’s hypothesis for the intercbangeabiliU' of space 
and time derivatives at low turbulence levels has been e.x- 
pressed in terms of criteria which depend upon turbulence 
level, Reynolds number, and X/X,. Applied to the flows 
studied here it shows that in these cases such a transforma- 
tion is permissible. By substitution of the theoretical e.x- 

pression for ^ (/?x), a slightly sim])ler and rougher criterion 

is (lerivod, depouding only upon turbulence level and /?x. 

The floHNs Hopkins University, 

Baltimore, Md., June 5, 1951, 
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TABLE I 

RESULTS FOR LAGRANGIAN AND EULERIAN 
SCALES AND MICROSCALES 
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1 
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8. 5 

1 

43.4 
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0.41 

0. 36 

0.28 

0. 43 
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1.48 
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H 

43.4 

2.0 
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.245 

.15 

.34 

1.48 

2.3 

43.5 

40 

.0136 


V 4 

43.4 

2.0 

.12 

. 12 

.08 

.21 

1.00 

2.6 

32 

21 

.0202 

25.6 

1 

86. 1 

1.5 

0. 34 

0. 22 

0. 39 

0.23 

0. 65 

0.59 

67 

77 
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86 1 

1.4 

.23 

.235 

.20 

.25 

1.02 

1.25 

42.5 

37 

.0138 


H 

86.1 

1.3 

.17 

. 11 

.10 

.19 

.65 

1.9 

29 

17 

.0203 

25.6 

H 

43.4 
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1.48 

2.3 

43.5 

40 
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.23 

.235 

.20 
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1.25 

42.5 

37 
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.22 

.30 
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49 
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.71 
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19 
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